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Abstract 

A general proecological urban road traffic control idea for the street canyon is pro- 
posed with emphasis placed on development of advanced continuum field gasdynamical 
(hydrodynamical) control model of the street canyon. The continuum field model of 
optimal control of street canyon is studied. The mathematical physics' approach (Eu- 
lerian approach) to vehicular movement, to pollutants' emission, and to pollutants' 
dynamics is used. The rigorous mathematical model is presented, using gasdynamical 
(hydrodynamical) theory for both air constituents and vehicles, including many types 
of vehicles and many types of pollutant (exhaust gases) emitted from vehicles. The six 
optimal control problems are formulated. 

1 Description of the model. 

In the present article we develop a continuum field model of the street canyon. In the 
next article we will deal with numerical examples [1]. The vehicular flow in the canyon 
is multilane bidirectional one-level rectilinear, and it is considered with two coordinated 
signalized junctions [2, 3, 4]. The vehicles belong to different vehicular classes: passenger 
cars, and trucks. Emissions from the vehicles are based on technical measurements and 
many types of pollutants are considered (carbon monoxide CO, hydrocarbons HC, nitrogen 
oxides NO x ). The vehicular dynamics is based on a hydrodynamical approach [5]. The 
governing equations are the continuity equation for the number of vehicles, and Greenshields' 
equilibrium speed-density u-k model [6]. 

The model of dynamics of pollutants is also hydrodynamical. The model consists of a set 
of mutually interconnected nonlinear, three-dimensional, time-dependent, partial differential 
equations with nonzero right-hand sides (sources), and of boundary and of initial problem. 
The pollutants, oxygen, and the remaining gaseous constituents of air, are treated as mixture 
of noninteracting, Newtonian, viscous fluid (perfect or ideal gases). The complete model 
incorporates as variables the following fields: density of the mixture, mass concentrations of 
constituents of the mixture, velocity of mixture, temperature of mixture, pressure of mixture, 
intrinsic (internal) energy of mixture, densities of vehicles, and velocities of vehicles. The 
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model is based on the assumption of local laws of balance (conservation) of: mass of the 
mixture, masses of its constituents, momentum and energy of the mixture, the numbers of 
the vehicles, as well as of the state equations (Clapeyron's law and Greenshields' model). 
The equations of dynamics are solved by the finite difference scheme. 

The six separate monocriterial optimization problems are formulated by defining the 
functionals of total travel time, of global emissions of pollutants, and of global concentrations 
of pollutants, both in the studied street canyon and in its two nearest neighbour substitute 
canyons. The vector of control is a five-tuple composed of two cycle times, two green times, 
and one offset time between the traffic lights. The optimal control problem consists of 
minimization of the six functionals over the admissible control domain. 



2 Assumptions. 

The geometrical assumptions of the model are as follows: 
Gl. The street canyon is represented by a cuboid: 

Q = [0,a] x [0,6] x [0,c], (2.1) 

with the boundary dQ composed of six walls: 

VI 

dQ = (J dflj. (2.2) 
j=i 

The structure of the canyon is simplified by the assumption that the walls of the buildings 
and the road surface are rectangles. If we put the origin of the Cartesian coordinate system 
(x,y,z) in cuboid's corner, then we have two canyon walls f2rv,fiv, at y — 0, y — b, and 
road's surface flyi, at z — 0. The remaining three non-solid open surfaces of air Qi, f2n, f^m, 
have the coordinates x = 0, x = a, z = c, respectively. We assume that always: 

(x,y,z)en, (2.3) 

and 

(2.4) 



dfli = 


= {(x,y,z) G Q : 


x = 


= 0}, 


dn u : 


= {(x,y,z) G Q : 


X ■ 


= a}, 


dtt m 


= {(x,y,z) G Q 


: z 


= c}, 


dQ w 


= {(x,y,z) G Q 


■ y 


= 0}, 


on y : 


= {(x,y,z) G Q : 


y - 


= b}, 




= {(x,y,z) G Q 


: z 


= 0}. 



G2. There are neither holes in the walls nor vegetation alongside the road. The remain- 
ing three surfaces of the cuboid also do not have holes since they simulate non-solid open 
rectangles of air. 



2 



G3. The road sections which constitute the bottom of the street canyon are rectilinear. 

G4. At each end of the street canyon there are entrance and exit junctions (M = 2) 
with traffic signals (their coordinates are x — 0, x — a). 

G5. The vehicles of VT distinguishable emission types are material points. The vehicles 
are treated as hydrodynamical fluid. There are ul = n\ left lanes and n 2 = ur right ones 
(the traffic is bidirectional). 

The physical assumptions of this model are as follows (cf. papers [7] - [50]): 

PI. The considered mixture of gases consists of N = Ne — 1 + Na gases. The first 
Ne — 1 = 3 gases are the exhaust gases emitted by vehicle engines during combustion 
(CO,CH,NO x , we neglect the presence of S0 2 ). The remaining Na = 9 gases are the 
constituents of air: 2 , N 2 , Ar, C0 2 , Ne, He, Kr, Xe, H 2 (we neglect the presence of H 2 0, 3 ). 

P2. The walls of the canyon and the surface of the road are impervious for all gases of 
the mixture. The remaining three surfaces of the cuboid are pervious for external fluxes of 
exhaust gases and air constituents. 

P3. The internal sources of air constituents are not present with the exception of oxygen, 
i.e., N E constituent of the gaseous mixture. There are internal mobile sources of exhaust 
gases (passenger cars and trucks, with many types of engines: diesel or petrol, and with 
different ages of engines). During combustion, the engine consumes oxygen, therefore with 
each internal mobile source of exhaust gases, a negative source of oxygen (sink) is connected. 

P4. The gaseous mixture is treated as a compressible, Newtonian, and viscous fluid. We 
assume that also the constituents of the mixture are compressible, Newtonian, and viscous 
fluids. The constituents do not interact with each other. The ith constituent possesses 
individual velocity Vj(x, y, z, t), density pi(x,y,z,t), and pressure Pi(x,y,z,t), whereas the 
mixture possesses total velocity v(x, y, z, t), densityp(x, y, z, t), and pressure p(x, y, z, t). We 
assume that 

v(x, y, z, t) = 2^ — , • Vi(x,y,z,t), 2.5 

~{ p(x,y,z,t) 

N 

p(x, y, z,t) = J2 Pi( x > V, z i *)> ( 2 - 6 ) 
i=i 

N 

P(x,y,z,t) = ^2pi(x,y,z,t). (2.7) 
i=i 

In order to simplify the set of equations governing the dynamics of mixture, we assume that 
the total velocity is equal to the velocities of the constituents 

v(x, y, z, t) = Vi(x, y,z,t),i = l,..., N. (2.8) 

Hence, we can restrict our attention to the equations of balance of total momentum of 
mixture El, of the total mass of mixture E2, of masses of constituents E3, and of the energy 
of mixture E4. We assume that the equation of state for mixture E5 is averaged over the 
constituents. We also consider equations of balances of the numbers of vehicles E6, as well 
as equations of state for vehicles E7. 



3 



3 Variables. 



The following set of descriptive dynamic model variables A0-A10 together with their bound- 
ary B0-B8 (for t > 0), and initial conditions C0-C7 (for t — 0), and with the set of equations 
E1-E8 that governs their dynamics, is assumed [2]. We always consider 

(x 1 j /l 2,t)es,s = flx[o 1 r s ] ) (3.1) 

where £ is manifold of domains of the fields, T s > is time of simulation. The border <9E 
of £ is composed of six subsets <9£j: 

VI 

<9£ = (J <9Ej, (3.2) 
j=i 

dEj = dQjX [0,T S ], J = I,...,VI. 

AO. T(x,y, z,t), temperature of the gaseous mixture. 

Al. v(x,y, z,t), total velocity of the gaseous mixture (compare Eq. (2.5)). 

A2. p(x,y,z,t), total density of the gaseous mixture (compare Eq. (2.6)). 

A3. Ci(x,y, z,t), mass concentration of the ith constituent of gaseous mixture, 

c i (x,y,z,t) = *j X ' V ' Z 'tl i = l,...,N (3.3) 
p(x,y,z,t) 



(let us note that due to the condition 



N 

^2ci(x,y,z,t) = l, (3.4) 



i=l 



one concentration of the constituent is a dependent variable). 

A4a. p(x,y, z,t), pressure of gaseous mixture (compare Eqs (2.7), (7.16)). 

A4b. pi(x,y, z,t),i = 1, N, partial pressure of the ith constituent of gaseous mixture 
(compare Eqs. (2.7), (7.16), (7.17)). 

A5. kf vt (x,t), density of vehicles of type vt on the Ith lane measured in [— ], where 
for n\ = riL lanes s — 1 and I is the left lane's number, I = 1, ...,riL, whereas for 112 = tir 
lanes s = 2 and / is the right lane's number, I = 1, ...,tir, vt is the vehicular type number, 
vt = 1,...,VT. 

A6. wf vt (x,t), velocity of vehicles of type vt on the Ith lane. 

A7. ef ctvt (x,t), emissivity of ctth constituent of exhaust gases from vehicles of type vt 
on Ith lane measured in [— 1, ct is number of constituent, ct — 1, CT. 

A8. u m = (g m ,C m , F), vector of control on the mth junction, m = 1,...,M (M = 2), 
which contains traffic signals green times g m , and cycle times C m , and offset time F between 
the traffic signals. The vector of control u reads: 

VL=(g 1 ,C 1 ,g 2 ,C 2 ,F). (3.5) 
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The admissible control domain set £/ adm for this vector in the simulation time period Th- 
reads: 

U* ba = {(g 1 ,C 1 ,g 2 ,C 2 ,F i ): (3.6) 

Qrn G i^9m,mu\i Qm,max) i Cm £ {C m ^ m \ a ., Cm,m&x) j -F £ (F m ; n , ^max) j 1; i 

whereas g m ,max — C m — g mt orth, where g mi0 rth are green times on orthogonal direction of the 
junctions (on the canyons orthogonal to the one studied), F max = C 2 — 5p, where 5p is unit 
step in direction of F in parameter space (compare FO). 

A9. G™ ut (g m , C m , F, t), traffic signal on mth junction at x — x m . For x\ — a the traffic 
signal G\ ut governs all left lanes (outgoing vehicles) and all right lanes (incoming vehicles), 
whereas for x 2 = the traffic signal G^ ut governs all right lanes (outgoing vehicles) and all 
left lanes (incoming vehicles). For the signals we assume the Boolean output values: GREEN 
and RED. 

A10. of t (x, t), rate of change of linear density of energy connected to heat produced by 
engines of vehicles of type vt on Ith lane measured in [^] . 



4 Boundary conditions. 

We assume the following boundary conditions [2]: 
BOa-BOc. T| aSx =T K \ T > K - 
BOd-BOf. V n T| 9SL = 0. 
Bla-Blc. v| 9Sk = v x | s ^. 
Bld-Blf. V n v| 9SL = K 
B2a-B2c. p\ dTiK = p K \x' K - 

B2d-B2f. V n p| 9Si =0. K 
B3a-B3c. = q K L> . 

B3d-B3f. V n Q|a SL = 0. 
B4a-B4c. p\ dT , K = p K \^. 

B4d-B4f. V nP |as L =0. K 
B5a-B5d. kf >vt (x m ,t) = kl vt , P (t). 
B6a-B6d. wf M (x m ,t)=wf >vt>P (t). 
B7a-B7d. el ct<vt (x m ,t) =e s lMP (t). 
B8a-B8d. af tVt {x m ,t) = <rf tVttP {t). 
We define additional sets: 

Yt I = {(y,z,t):{x,y,z,t)ei:} = Y! Il , (4.1) 
s m = {(x,y,t) : (x,y,z,t) G £}, 

and we assume that K — I, II, III, L = IV, V, VI. The gradient operator V n works in direction 
of unit normal outward vector n to border <9S^. V n v is gradient of vector (so it is a tensor of 
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rank 2), O is a zero tensor. P = in, out, is the input and output index, and we have the fol- 
lowing combinations of triads of indices: (s, P, m) = (1, in, 2), (2, in, 1), (1, out, 1), (2, out, 2), 
respectively for B5-B8. Conditions Bld-Blf result from viscosity of the gaseous mixture 
since the velocity of viscous fluid on immobile and impervious surface is zero. Similarly, con- 
ditions BOd-BOf, B2d-B2f, B3d-B3f, B4d-B4f result from the fact that the walls and the 
surface of the road are impervious solid bodies. According to [5], we assume the boundary 
conditions B5a-B5d in the form: 
B5aS-B5dS. 

Kvt,p(t) = ^arrival, if GZt{C m , 9m, F, t) = GREEN, and QUEUED) = FALSE, 
KvtA*) = *W, if GZt(C m ,9m,F,t) = GREEN, and QUEUE(x m ) = TRUE, 
k>t,vt,p{t) — ^%t,jam; if G™ ut (C m , g rn , F, t) = RED, 

where QUEUE(x m ) = TRUE/FALSE means that there exist/does not exist a queue 
at x = x m , and ^„ tarrival , kf^ t , and k s lvt - J&m , are arrival, saturation, and jam vehicular 
densities, respectively (compare Tables 3, and 8 of [1]). 

The existence of the queues at the entrances to the canyon (at x 2 = for the left lanes, 
and at X\ — a for the right lanes) is determined by the values of the vehicular densities 
changing in the following way: 

B5aSS-B5dSS. 



kt vt (Z s ,t) = kl vt>Q MteA s Q , (4.2) 

+oo 

^GREEN = 

[0,T s ]n (J [* s + nCi,* a + nCi + ^i), (4.3) 

n=— oo 

A S RED = [0, T S ] - ^ REEN , Q = GREEN, RED, (4.4) 
Zi = -S x ,Z2 = a + 6 x ,t 1 = 0,t 2 = F, (4.5) 



where fcf„ t green > k^ vtRED , are green and red vehicular densities, respectively [1], and 5 X is 
the unit step in x-direction in the domain space E. 

Remark: The functions: Tj, vj, pj, Ci : j, pj, kf vt P , wf^p, e s l ct vt P , <jf vt P , are given and 
they fulfill the natural constraints: 

N 

J2c iyJ (x,y,z,t) = 1. (4.6) 

i=l 

5 Initial conditions. 

The initial conditions are as follows [2]: 
CO. T|s = Tq\q. 
CI. v|s = v |n- 
C2. p| So = pob- 
C3. q|s = Cj ; o|n- 
C4. p| So =po|n- 
C5a-C5b. kl vt (x,0) = kf tVtfi (x). 
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C6a-C6b. wl vt (x,0) = wl vtfi (x). 
C7a-C7b. el ctM (x,0) = ef MO (x). 
C8a-C8b. af tVt (x,0) = <rf tVt>0 (x). 

Remark: The functions: T , v , po, c ij0 ,Po, tf M ,o, w f,«t,o> e ? )C t,«t,o> ^%t,o> are g iven and the y 
fulfill the natural constraint: 

£ci,o(a:,2/,2,f) = 1. (5.1) 
i=i 

6 Sources. 

In order to represent the emission process, we assume the following internal sources [2]: 

DO. a(x,y, z,t), the rate of change of the volume density of internal sources of energy 
connected with the production of heat by vehicular engines, measured in [^-] . We assume 
that the sources of energy are situated in n s left and right lanes at y — yf, at the level of the 
road z = 0: 

I 2 n s VT 

a(x,y,z,t) = - 2 a tvt(x^)XD^{x,y,z), (6.1) 

' C s=l 1=1 vt=\ 

where 

Dt = {(x,y,0) : (x,y,0) etl,y = yf}, (6.2) 
are the vehicular lanes, and 

/ \ f 1 ior(x,y,z)eD ,„ 

x^y^) = { f r(,;L)^ ' (6 - 3) 

is the characteristic function of set D. 

Dl. S(x,y, z,t), the rate of change of the volume density of internal sources of gaseous 
mixture consisting of exhaust gases and oxygen, measured in [— ]. 

D2. S^ t (x,y, z,t), the rate of change of the volume density of internal sources (the 
emission rate) of the ctth constituent of exhaust gases emitted by all vehicles in the canyon, 
measured in [^]. We assume that the sources of exhaust gases are situated in n s left and 
right lanes y — yf, at the level of the road z = 0: 

Y 2 n 3 VT 

S%(x,y,z,t) = r--EEE e i,ct,vt( x ^)^(x,y,z). (6.4) 

6 s=l 1=1 vt=l 

Sh (x,y, z,t), the volume density of negative internal sources (the emission rate) of 
oxygen absorbed by all vehicles in the canyon, measured in [^]. We assume that 

S§ E (x, y, z, t) = ONOX • S^ix, y, z, t), 
where ONOX = -0.5308. The following relation holds: 

N E 

S(x,y,z,t) = J2 S%,(x,y,z,t). (6.5) 

ne=l 
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7 Equations of dynamics. 



Under the above model specifications, the complete set of equations of dynamics of the model 
is formulated as follows (we follow the general idea presented in [7, 8]): 

El. Balance of momentum of mixture - Navier Stokes equation. 

P(~ft + ( v ° V ) v ) + Sv = -^P + »?Av + (f + ^)V(divv) + F, (7.1) 

where 77 is the first viscosity coefficient (77 = 18.1 • 10~ 6 [^] for air at temperature T = 
293.16 [K]), £ is the second viscosity coefficient (£ = 15.6 • 10 _6 [^] for air at temperature 
T = 293.16[K], compare [51]), F = pg is the gravitational body force density, g is the 
gravitational acceleration of Earth (g = (0, 0, — 9.81)[jj|]), Vv is gradient of the vector (so 
it is a tensor of rank 2). We assume that the gaseous mixture is a compressible and viscous 
fluid. 

E2. Balance of mass of mixture - Equation of continuity. 

^ + div(pv) = S. (7.2) 

We have assumed the source Dl. 

E3. Balances of masses of constituents of mixture - Diffusion equations. 
E3a. 

p (-l +V oVc t ) = S?-c t S+ (7.3) 

N-l u 

+ E {(Am - D iN ) ■ div[pV(c m + -^VT)]}, i = 1, .., N E . 



T 

m=l ± 



E3b. 



p(^ + voVq) = -c^+ (7.4) 

N-l u 

+ E {(Am - D iN ) ■ div[ P V(c m + ^VT)]}, i — (N E + 1), N, 

m=l 1 

where Di m = D m i is the mutual diffusivity coefficient from the i-th constituent to m-th one, 
and Da is the autodiffusivity coefficient of the i-th constituent, and kr^m is the thermodif- 
fusion ratio of the m-th constituent. The diffusivity coefficients and thermodiffusion ratios 
are constant and known (compare [52]). In E3a we have assumed the sources D1-D2. In 
E3b only the source Dl is taken into account. Since the mixture is in motion, we cannot 
neglect the convection term: vo Vq. We assume that the barodiffusion and gravitodiffusion 
coefficients are equal to zero. 

E4. Balance of energy of mixture. 

de 1 
P<K dt + v Ve ) = ~(~2 ^ + e ) 5 + T : Vv + div(-q) + (7, (7.5) 
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where e is the mass density of intrinsic (internal) energy of the air mixture, T is the stress 
tensor, symbol : denotes the contraction operation, q is the vector of flux of heat. We assume 
that [2]: 

e = f>, ( 7 - 6 ) 

i=i 

e t = _ {ci fc fl Texp(-— ) • [(--) • (1 - exp(-— )) - exp(-— )]} +(7.7) 

+/^iCjj 

fii = — , (7-8) 

mi 



c p,i 777 ■ O.Tfh ^ 

Hi = k B T ■ {\n[(c iP )(k B T)-^(^)( )*]} + mi \g\z, (7.9) 



air \ / — />»( 1, 

mi ' mi 

T mfc = - P 5 mfc + (7.10) 
+17 • [(|^ + - ^ d ^ v ))|rH + * ' [(^div(v)) 2 ], m, fc = 1, 3, 

T:Vv=EET mfc |^, (7.11) 

m=l fe=l 

q = E{[(^ + ^] + [(-«) VT]}, (7.12) 

1=1 "»» 

j l = - P A,(Q + ^VT), (7.13) 
a .. = (7 14) 

A = [pA,] • { k -f - tM^^ -A , (7.15) 

i\ dCi '( c n)n = l,...,N,i^n,T,p\ 

where q is the mass density of intrinsic (internal) energy of the ith constituent of the air 
mixture, rrii the molecular mass of the ith constituent, k B = 1.3807- 10~ 23 [^] is Boltzmann's 
constant, /z, is the complete partial chemical potential of the ith constituent of the air mixture 
(it is complete since it is composed of chemical potential without external force field and of 
external potential), m a i r = 28.966 [u] is the molecular mass of air (l[u] = 1.66054- 10 -27 [kg]), 
5 m k is Kronecker's delta, c Pt i is the specific heat at constant pressure of the ith constituent of 
air mixture, h = 6.62608 • 10~ 34 [J ■ s] is Planck's constant, jj is the vector of flux of mass of 
the ith constituent of the air mixture, and k is the coefficient of thermal conductivity of air. 
These magnitudes are derived from Grand Canonical ensemble with external gravitational 
Newtonian field. 

E5. Equation of state of the mixture - Constitutive equation - Clapeyron's 
equation. 



9 



P rn air 

is Clapeyron's equation of state for a gaseous mixture, where R = 8.3145 [ mo / e . K ] is the gas 
constant. 

Pi = Ci-— BE --p (7.17) 

rrii 

are partial pressures of constituents according to Dalton's law. 

E6. Balances of numbers of vehicles - Equations of continuity of vehicles. 



dkf vt 



+ dw(kf vt wf vt ) = 0. (7.18) 



dt 

E7. Equations of state of vehicles - Greenshields model. 

wUM) = • (! - %^),0,0). (7.19) 

Z,ui,jam 

The Greenshields equilibrium speed-density u-k model is assumed [6]. The values of 
maximum free flow speed wf vt j, and of jam vehicular densities kf vt ^ am , are given in Tables 
3 and 8 of [1]. 

E8. Technical parameters. 

The dependence of emissivity on the density and velocity of vehicles is assumed in the 
form [54]: 
E8a. 



ei,ct,vt( x > *) = fytrffo *) ' K-r ^ ) • {e c t,vt,i l+ i - ectMA) + e<*,vt,ii\, ( 7 -20) 

w ct,vt,ii+l ~ w ct,vt,ii 

where w rfi „ t)i; are experimental velocities, \wf M (x,t)\ G {wct,vt,iv Wct^u+i), &ct,vt,i v are ex- 
perimental emissions of the ctth exhaust gas from single vehicle of vtth type at velocities 
Wct iV t,i v respectively, measured in [-^-], %i = 1, Nem, Nem is the number of experimental 
measurements. Similarly, the dependence of the change of the linear density of energy on 
the density and velocity of vehicles is taken in the form: 
E8b. 

0j%t(z> t) = q vt ■ kl vt (x, t) ■ [( Wl ! vt X,t _ Wvt ' H ) . (a vt:il+1 - a vtM ) + a vtjil ], (7.21) 

w vt,ii+l ~ w vt,i t 

where a vt ^ n are experimental values of consumption of gasoline/diesel for a single vehicle of 
vtth type at velocities w vt:iv respectively, measured in [-^-], q v t is the emitted combustion 
energy per unit mass of gasoline/diesel [^] (compare [53]). 
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8 Optimization problems. 

Our control task is the minimization of the measures of the total travel time (TTT) [5], 
emissions (E), and concentrations (C) of exhaust gases in the street canyon, therefore the 
appropriate optimization problems may be formulated as follows [2]: 
FO. Vector of control. 

u=(g 1 ,C 1 ,g 2 ,C 2 ,F)eU adm , (8.1) 

where u is vector of boundary control, g m are green times, C m are cycle times, F is offset 
time, and JJ adm is a set of admissible control variables (compare A8, A9, B5, B5S, B5SS). 

We define six functional F1-F6 of the total travel time, emissions, and concentrations 
of pollutants in single canyon, and in canyon with the nearest neighbour substitute canyons, 
respectively. 

Fl. Total travel time for a single canyon. 

2 n s VT a Tg 



* '"s v pa pis 

Jttt(u) = EEE/ / KvtMdxdt. (8.2) 

s=l 1=1 vt=l J0 J0 



F2. Global emission for a single canyon. 

2 n s CT VT a Tg 

^e(u)=EEEE/ / el ct Jx,t)dxdt. (8.3) 

s=l 1=1 ct=l vt=l J ° J ° 

F3. Global pollutants concentration for a single canyon. 

J c (u)=pstp- E / / / / Ci(x,y,z,t)dxdydzdt. (8.4) 
~[ Jo Jo Jo Jo 

F4. Total travel time for the canyon in street subnetwork. 

^TTT,cxt(u) = JtTt(u) + (8.5) 

2 n s VT 

+a ■ ^2 a TTT,ext E E Kvtjam ' fts ~ 9s)- 
s=l 1=1 vt=l 

F5. Global emission for the canyon in street subnetwork. 

J E) ext(u) = J E (u) + (8.6) 

2 n s CT VT 

+ a ■ E «E,cxt E E E <rf,^,jam '(C 8 -g s ). 
s=l 1=1 ct=l vt=l 

F6. Global pollutants concentration for the canyon in street subnetwork. 

J C ,ext(u) = J C (U) + (8.7) 

N E -1 2 

+Pstp -a-b-c- c ^stp • E a c,cxt • (C s - 9s)- 

i=l s=l 
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The integrands kf vt , ef ct vt , ci in functionals F1-F6 depend on the control vector u FO 
through the boundary conditions B0-B8, through the equations of dynamics E1-E8, as 
well as, through the sources D0-D2. The value of the vector of control u directly affects 
the boundary conditions B5, B5S, B5SS, and then the boundary conditions B6-B8 for 
vehicular densities, velocities, and emissivities. It also affects the sources D0-D2. Next, 
it propagates to the equations of dynamics E1-E8 and then it influences the values of 
functionals F1-F6. We only deal with six monocriterial optimization problems 01-06, 
and not with one multicriterial problem. We put the scaling parameters equal to unity: 
a TTT,ext = a E,ext = a c,ext = 1-0, in functionals F4-F6. pstp is the density of air at standard 
temperature and pressure STP, q 5 stp is concentration of the ith constituent of air at standard 
temperature and pressure. Jttt and JrTT.ext are measured in [veh • s], J E and J EiCx t are 
measured in [kg], and Jq and Jc,ext are measured in [kg ■ s], respectively. 

Now we formulate six separate monocriterial optimization problems 01-06 that consist 
in minimization of functionals F1-F6 with respect to control vector FO over admissible 
domain, while the equations of dynamics E1-E8 are fulfilled. 

01. Minimization of total travel time for a single canyon. 

J* TT = J TTT (u^ TT ) = min{u G C/ adm : J T tt(u)}; (8.8) 

02. Minimization of global emission for a single canyon. 

J E = J E (u E ) = min{u G f/ adm : J E (u)}; (8.9) 

03. Minimization of global pollutants concentration for a single canyon. 

r c = JcK) = min{u G C/ adm : J c (u)}; (8.10) 

04. Minimization of total travel time for a canyon in street subnetwork. 

4rr,«t = ^TTT,cxt(uTTT,cxt) = min{u G f/ adm : J TTT ,cxt(u)}; (8.11) 

05. Minimization of global emission for a canyon in street subnetwork. 

^E,ext = </E,cxt« cxt ) = min{u G C/ adm : J E , cxt (u)}; (8.12) 

06. Minimization of global pollutants concentration for a canyon in street 
subnetwork. 

^,cxt = </c,cxt« ext ) = min{u G f/ adm : J c ,ext(u)}, (8.13) 
where J^tt) ^ej ^hi "^ttt ext 1 ^Eext) ^Cext; are the minimal values of the functionals F1-F6, 
and u^ TT , u E , u^, u^TT.ext) u E,ext5 u c,cxt5 are control vectors at which the functionals reach 
the minima, respectively. 
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